Fermionic superfluidity with positive scattering length 
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Superfluidity in an ultracold Fermi gas is usually associated with either a negative scattering 
length, or the presence of a two-body bound state. We show that none of these ingredients is neces- 
sary to achieve superfluidity. Using a narrow Feshbach resonance with strong repulsive background 
interactions, the effective interactions can be repulsive for small energies and attractive for ener- 
gies around the Fermi energy, similar to the effective interactions between electrons in a metallic 
superconductor. This can result in BCS-type superfluidity while the scattering length is positive. 

PACS numbers: 03.75.Ss, 34.50.-s, 64.60.-i, 74.20.-z 
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The crossover between Bardccn-Coopcr-Schricffcr type 
superfluidity (BCS) and Bose-Einstein condensation 
(BEC) is a fundamental and challenging problem in 
physics. In ultracold Fermi gases, both types of super- 
fluidity are related to the condensation of atom pairs, 
which can be regarded as composite bosons. A crossover 
between the BCS and BEC regimes, which has beenpre- 
dicted theoretically already several decades ago P, 0, d , 
can be achieved by changing the strength of the inter- 
actions. While several experiments have studied certain 
aspects of this strongly interacting state 0, IE IE j 0Iu y 
recently the existence of superfluidity in the crossover 
regime has been unambiguously demonstrated by exper- 
iment @. 

Experimentally the crossover is realized by changing 
the strength of the two-body interactions via an s-wave 
Feshbach resonance |3, |lfj . A measure for this interac- 
tion strength at zero collision energy is provided by the 
scattering length a, which is the first parameter in the ef- 
fective range expansion of the scattering phase shift • 
It has been suggested that the physical properties in the 
crossover only depend on the parameter (fcpa) -1 , with 
kp the Fermi wavenumber. This has the strong impli- 
cation that exactly on resonance, where a goes through 
infinity, all thermodynamic and superfluid properties are 
independent of the microscopic details of the interactions, 
and therefore should be universal quantities for all atomic 
systems under study p^3. fl3L ll4L fl5|| . 

However, the microscopic physics that gives rise to the 
Feshbach resonance reveals that more than one parame- 
ter is needed to describe the scattering phase shift. Apart 
from the scattering length, the width of the Feshbach 
resonance is the most important parameter. In order to 
understand the origin of this width, we have to sepa- 
rate the resonance from the background interactions. At 
large relative distance, the two interacting atoms are sen- 
sitive only to the background interaction potential. At 
short distance, this potential is coupled to another poten- 
tial with a different internal spin configuration. The lat- 
ter potential is energetically closed at long distance, and 
exhibits a molecular state close to the collision thresh- 
old |E M- 

The coupling strength between these two 
potentials can be related to an energy width which, in 
comparison to the Fermi energy of the system, catego- 



rizes a Feshbach resonance in terms of its broadness. This 
width can be converted into an effective range coefficient, 
which is the second parameter in the effective range ex- 
pansion ^lj. Using this approach, it has been pointed 
out in several papers that universality does not hold for 
narrow Feshbach resonances 0, 0, 0, . 

It should be noted that in these papers the background 
interactions were taken to be small, and the Feshbach 
resonance was described using an effective range coeffi- 
cient independent of magnetic field or momentum. How- 
ever, in case the background interactions do have a reso- 
nant character - which is indicated by a large background 
scattering length, as for example in some resonances in 
6 Li, 133 Cs, and 85 Rb - the scattering phase shift shows 
non-trivial energy dependence beyond the effective range 
approximation, due to the inter play between the back- 
ground and Feshbach resonance |2Cj . 

In this paper we investigate the BCS-BEC crossover 
for a Feshbach resonance which allows for a large back- 
ground scattering length. We show that we can mimic 
a situation well known from condensed-matter physics, 
where Cooper pairs can be formed due to an effectively 
attractive interaction for relative wavenumbers around 
kp, while for smaller wavenumbers the interaction is ef- 
fectively repulsive. For electrons in a metal such an inter- 
action exists due to the combined effect of the Coulomb 
and electron-phonon interactions plj . while in atomic 
systems this can be achieved by the interplay of a Fesh- 
bach and a background resonance. 



I. FESHBACH RESONANCES WITH LARGE 
BACKGROUND SCATTERING LENGTH 

Feshbach resonances are presently an indispensable 
tool to control the atomic interactions in cold atomic 
physics. By simply changing the magnetic field the in- 
teractions can be tuned from weak to strong, and from 
effectively attractive to repulsive. We consider a two- 
component Fermi gas in the dilute limit, i.e. the density 
n satisfies nr^ <C 1, where r ~ (mCe/ 16H 2 ) 1 / 4 [^E] is 
the range of the potential, with m the atomic mass and 
C§ the van der Waals coefficient. The Fermi wavenumber 
kp is related to the density of a single spin component 
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as n j = ni = n — kp/6iT 2 . 

Feshbach resonances originate from the coupling of 
scattering atoms in an energetically open channel, to one 
(or several) bound state(s) in energetically closed chan- 
nels. When the energy of a closed-channel bound state 
is close to the energy of the scattering atoms, the scat- 
tering process becomes resonant. The effect of two-body 
interactions on the atomic wave functions can be sum- 
marized by the energy-dependent scattering phase shift 
S(k), where k — ymE/F, with E the relative collision 
energy. The scattering matrix is related to the phase 
shift by S(k) — exp2j<5(fc), and for the cold alkali gases 
under consideration its most general form is given by pcj 



From Eq. J3J it follows that the total scattering length is 
related to these parameters by 



S(k) = Sp(k) 



E - ^we - A(£) - |TCg) 
E - t/ bare - A(E) + iT(E) ' 



(1) 



The direct part of the scattering matrix Sp(k) describes 
the interactions in the open channel, and A(E) — ^T(E) is 
the complex energy shift, which describes the dressing of 
the closed-channel bound state with bare detuning Vbme 
due to the coupling with the open-channel atoms. 

In scattering systems with a background scattering 
length Ob g — 7-o, the direct interactions are non-resonant, 
and the energy shift takes the simple form A — iCk. The 
bare detuning is dressed by the constant term A such 
that the dressed detuning is given by v = fbare + A. The 
scattering length diverges for v = 0. The term iCk de- 
scribes the finite lifetime of the bound state coupled to 
the continuum, where the coupling strength C is a con- 
stant related to the magnetic field width of the resonance 
(see below). 

In scattering systems with a background scattering 
length |&bg| ^> 7"o, the open channel in general has a 
bound or virtual (nearly-bound) state close to the zero- 
energy collision threshold. In this case, the resonant en- 
ergy dependence of the open-channel propagator gives 
rise to non-trivial energy dependence of the energy shift. 
Explicitly taking into account the resonance poles of both 
the closed-channel bound state and the open-channel 
bound (or virtual) state, the energy shift can be written 
as A(E) - |r(£) = A - iCkQ 1 HJ. Here the com- 
plex energy shift is written as its value at zero collision 
energy A = A(E = 0), which is real, and an additional 
term that accounts for the energy dependence of both 
the real part of the complex energy shift A(E) — A, as 
well as the imaginary part T(E). The coupling strength 
C is related to the magnetic field width AB = B' — Bq 



as C = A/z mag ABa i 



A/j ma «( J B - B ) is the dressed 



detuning, A/z mas is the magnetic moment difference be- 
tween the open and closed channel, and £fc = 1 + ika p 
with a p — abg — 7"o describes the additional energy de- 
pendence related to the resonant direct interactions. Bo 
is the field where a(B) diverges, and B' is the field where 
a(B) = ro- The interplay between the two resonances is 
thus described by the scattering phase shift: 

S(k) = -kr - argCfc - arg (E - v + iCkQ 1 ). (2) 



a(B) =r + , 



1 



AB 
B-B* 



(3) 



In the following, we will discuss scattering systems with a 
resonant direct interaction between the fermions. In this 
case the non-resonant interaction terms related to ro will 
generally have a small effect compared to the resonant 
terms, and will be ignored in the remainder of this paper. 

The system can be described equivalently by taking 
into account the scattering poles of the T matrix of 
both the molecular and background bound (or virtual) 
states 11], which can be shown to correspond to solu- 
tions of 

1/a + ik - R e (k,B)k 2 = 0, (4) 

where R e (k,B) = h 2 {l + ika p )/(ma p Afi nia z(B-B' )) is 
a generalized energy and field dependent effective range 
parameter. These poles determine the two-body T ma- 
trix and scatterin g ph ase shift uniquely. In contrast to 
previous work [Ta Il7l ITsl Il9| , where R e was taken to be 
a constant, here R e depends on both energy and magnetic 
field and describes the most general Feshbach resonance 
that can be found in ultracold alkali gases. Note that for 
magnetic fields where a ~ ro, the effective range param- 
eter R e actually diverges and cannot be neglected even 
in the case of a broad resonance. The usually considered 
expressions for Feshbach resonance scattering without an 
open-channel resonance |l(j can be seen as a limit of the 
general expression Eq. J2J) for small <2b g . 



II. MANY-BODY DESCRIPTION 

In order to describe the many-body properties of the 
system, we use the thermodynamic approach introduced 
by Nozieres and Schmitt-Rink (NSR) |3J . This approach 
is based on a simplified description of the full crossover 
problem and cannot be expected to give quantitatively 
precise results. Therefore, in recent years several direc- 
tions have been pursued in order to formulate a more 
complete crossover theory pll2l I2I I2II2III2I |29|. 
However, the NSR description does capture the essential 
physics giving rise to the crossover phenomenon. Es- 
pecially the very simple relation between the two-body 
phase shifts in vacuum, and the in-medium phase shift 
used in NSR make this approach a very convenient frame- 
work to study the role of energy dependent scattering in 
a many-body system. In the following, we show how to 
generalize this approach in order to include the energy- 
dependent phase-shift given by Eq. J2J). 

Within the NSR approach, the two-body interactions 
are taken into account in the ladder approximation. The 
two-body interactions are described by the many-body 
T matrix, which can be found by solving the Lippmann- 
Schwinger equation (LS) 



T = U - U X T, 



(5) 
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where U is the two-body interaction potential, that gen- 
erally depends on the relative collision energy E, and x 
is the two- particle propagator. In a vacuum, the many- 
body T matrix should reduce to the two-body T matrix, 
which is known analytically and gives rise to the scatter- 
ing phase shift given by Eq. (J2J). The bound states and 
resonances are described by the poles of the two-body 
T matrix, and uniquely characterize the scattering prop- 
erties in the two-body sector. We are thus faced with 
constructing a potential U that gives the correct poles in 
the two-body limit. Using this potential ensures that the 
two-body physics is built correctly into the many-body 
theory. 

It can be shown from coupled-channels scattering 
theory that the effective interaction potential for the 
fermions in the open channel is given by 



Uw{E) = Vw 



E - Vba 



(6) 



where V is related to the direct interactions between 
fermions in the open channel, and the second term de- 
scribes the coupling of fermions to and from the closed- 
channel bound state [3(| • To simplify the equations, it is 
convenient to replace the microscopic finite-range inter- 
action potentials with contact potentials that correctly 
reproduce the scattering properties described by Eq. (0) 
in the two-body limit. However, since contact poten- 
tials do not explicitly depend on Fourier momentum, they 
introduce unphysical ultraviolet divergences which have 
to be regulated explicitly by renormalizing the potential 
and two-body propagator. The appropriate regulariza- 
tion procedure has been discussed in detail in Ref. |3l| . 
where it was shown that the T matrix is described by the 
renormalized LS equation 



T(E) = U cS (E) - U cS (E)xo(E+)T(E), 



(7) 



where the renormalized effective interaction is related to 
the microscopic parameters of section I by 



U eS (E) 



AttH 2 



C 



E-v 



(8) 



and the two-body pair-propagator in vacuum is explicitly 
regularized to account for the ultraviolet divergences: 



1 



2e k - E^ 



1 

2^ 



(9) 



with E + approaching the real axis from above, and 
€k = h 2 k 2 /2m. It can be shown that Eq. Q with the 
effective potential U c g(E) Eq. @ exactly reproduces the 
pole equation of Eq. |@J|. 

The many-body T matrix is obtained by replacing the 
regulated pair-propagator in vacuum by the regulated in- 
medium two-particle propagator, which is defined as 



x(q^ + )=E 



2efc — 2e„ 2e^ 



(10) 



where the relative energy of the atoms E — 2e p defines a 

relative momentum p = -^p (e + — \e q + 2/j) 1 ^ 2 , with the 
total energy of the scattering atoms (relative to the chem- 
ical potential) e + approaching the real axis from above, 
and ^e q the center-of-mass momentum of the scattering 
particles. The medium effects are described by the Fermi 
distribution functions 



i 



°P(t±-n) +i' 



(11) 



with (3 the inverse temperature and £± = e q / 2 ±k- The in- 
teractions are thus described by the many-body T matrix 



, -l 



f/ off (2e p )~ 1 + X (q,e + ). (12) 



This way of replacing the bare interactions by physical 
(renormalized) interactions is quite general, and closely 
related to a method due to Gorkov (33 , and used subse- 
quently in different forms (3^, . 

We now proceed to give a description of the many- 
body system in the NSR approach. The number den- 
sity N/V = 2n, with N the total number of atoms 
in a volume V, is found at a temperature T (in the 
normal state) by first calculating the thermodynamic 
grand potential Q within the ladder approximation, and 
then taking the derivative to the chemical potential fi: 
2n(fi,T) = — -*fj2y- The grand potential is then found to 
be 



n = fi + i yVe)0(q, e), 

•7T ' 



(13) 



where fio corresponds to a non-interacting gas of 
fermions, g(e) = (exp(/3e) — 1) _1 is the Bose distribu- 
tion function, and the phase of the many-body T matrix 
4> = (f> p + <jfi has a contribution from the resonant back- 
ground 



(q,e) = Imlog 1 



as well as from the Feshbach resonance 



(q,e) = Imlog \v - 2e p 



4irh 2 C 



II(q,e+) , (15) 



where II (q, e+) - X (q, e + ) (l + (q, ' m- 

eludes the dressing of the Feshbach molecules due to 
the resonant background interactions. Note that we ex- 
clude the contribution from a possible pole associated 
with the background interactions when a p > in the 
energy integral in Eq. (jl3|) . since we assume that the 
corresponding molecular state cannot be populated. The 
critical temperature is determined from the gap equa- 
tion T _1 (q = 0, e = 0) = 0, which relates the instability 
of the gas towards the formation of condensed pairs to a 
divergence in the many-body T matrix. Solving the num- 
ber and gap equations self-consistently gives the critical 
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FIG. 1: (Color Online) Critical temperature T c in the 
crossover region as a function of (fc^a) -1 in the case of broad 
and narrow Feshbach resonances with different background 
interaction strengths k F a p . Dashed line: broad resonance 
where (kpa) -1 is the only relevant parameter; solid, dash- 
dotted and dotted lines: narrow resonance where C = 1 and 
kF<i p = +0.5, 0, —0.5 respectively. 

temperature T c and chemical potential fi c at this tem- 
perature as a function of the interaction parameters. 

In the vacuum limit, where the two-particle propaga- 
tor reduces to imp/ATrh 2 and the energy of the scattering 
particles is e = 2e k + \z q — 2fi, the phase shift 0(q, e) only 
depends on the relative energy 2e k and is related to the 
scattering phase shift Eq. © by </>(q, 2e k + ^e q — 2\i) = 
—S(k) — n. Since k F r <C 1, we neglect the range of the 
background potential, and then the effective interaction 
Eq. (jSJ exactly reproduces the non-trivial energy depen- 
dence of Eq. ©. 

III. RESULTS AND DISCUSSION 

We now turn to the numerical results. We first consider 
the case of a broad Feshbach resonance, characterized by 
C = Ak F C/3Tr^ F > 1 or equivalently \k F R e (0, B )\ < 1. 
We solved the number and gap equations numerically 
with C = 100, and for different values of —1 < k F a p < 1. 
The results are not sensitive to the actual value of k F a p , 
and an example is shown in Fig. ^ (dashed line). This 
reproduces the well known results previously obtained 
by NSR and others HO^. In this case the only relevant 
parameter in the crossover region is the total scattering 
length Eq. and details of the interactions at higher 
energy scales do not enter the phase integral Eq. as 
they are cut off by the distribution functions. 

In the case of a narrow Feshbach resonance (C < 1), 
details of the interactions turn out to be essential in the 
crossover region, as well as in the BCS- and BEC-limits. 
In Fig. ^ the critical temperature in the crossover re- 
gion is shown for several narrow resonances (C = 1) with 
different values of the background interaction strength 
k F a p (solid, dash-dotted and dotted lines). In the BEC 



FIG. 2: (Color online) Critical temperature T c in the BCS 
limit, for small negative and positive scattering length k F a, 
in the case of a narrow Feshbach resonance with k F a p = +1 
(solid line), and of a broad resonance (dashed line). Inset: 
chemical potential fi c at T c for identical parameters. 



limit (right side of figure), the critical temperature ap- 
proaches the limiting value T c ~ 0.218Tp, and super- 
fluidity is dominated by the condensation of molecules. 
Compared to a broad resonance, T c is suppressed in 
the crossover region, while in the BCS-limit (left side of 
figure), the critical temperature remains a much larger 
fraction of the Fermi temperature T F compared to the 
broad resonance case. Somewhat surprisingly, T c be- 
comes larger in the BCS-limit if the background inter- 
actions are more repulsive. 

This leads us to the very interesting situation when the 
background scattering length is large and positive, while 
the Feshbach resonance is narrow. In Fig. [21 the criti- 
cal temperature for a narrow resonance (C — 1) with 
strong repulsive background interactions (k F a p = +1) is 
shown (solid line) [2f|. The cusp around T C /T F = 0.1 is 
not physical, but an artifact resulting from the NSR ap- 
proach [8fil l8^ | . If the system is initially on the BCS-side 
of the resonance {k F a ~ —0.5), and subsequently tuned 
away from resonance towards the point where a becomes 
close to zero, we find that the critical temperature can 
still be a large fraction of the Fermi temperature. This 
in contrast with a broad resonance (dashed line), where 
the superfluid state is lost before a changes sign. Go- 
ing through the point a = 0, the low-energy interactions 
become effectively repulsive as a > 0. However, the fre- 
quency dependence of the phase shift can still give rise 
to attractive interactions close to the Fermi energy. In 
Fig-Elthe phase (f>(k) is shown for the same resonance pa- 
rameters as in Fig. |21 and with the total scattering length 
equal to k F a — +0.3. A positive phase corresponds to 
effectively repulsive interactions, while a negative phase 
corresponds to attractive interactions. In the case of a 
broad resonance, we immediately see that for a > the 
interactions are repulsive at all relevant energies. For the 
narrow resonance with strong repulsive background inter- 
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FIG. 3: (Color online) Phase shift <j>(k) = 0(q, 2e fe + |e a - 
2/i) for fcfo = +0.3. The resonance parameters are identical 
to those of Fig. [5] Inset: scattering length as a function of 
magnetic field. The shaded area indicates the region where 
a > 0, but there are no molecules in the system. 



actions, although the interactions are repulsive for low 
collision energies, the atoms colliding at energies close to 
the Fermi energy feel an attractive interaction, which is 
the reason that the system can still become superfluid. 

From the inset of Fig. [21 we see that for kpa > 0.35 
the chemical potential becomes larger than /i p , reflecting 
the fact that superfluidity can be retained even when the 
total system acts as effectively repulsive. The repulsive 
behavior is indeed expected from the positive scattering 
length for the fermions. The fact that this system with 
repulsive interactions can still form a superfluid is a direct 
consequence of the energy dependence of the interactions 
for a narrow Feshbach resonance. 

We have done an independent check of our calculation 
by comparing this interesting regime to the regular BCS 
limit. Since we are in the regime where T c /Tp < 0.1 and 
/i ~ (ip, the BCS result can be expected to become valid 
again. We compared the functional dependence of T c /Tp 
versus hp a by using the well-known BCS formula 



Instead of using the scattering length a, we replace it by a 
'energy-dependent scattering length' a(p), related to the 
relative collision energy via p — WrnE/W, defined from 
the scattering phase shift Eq. (1): 



a(p) = -tan[f5(p)] /p. 



(17) 



T C /T F = - exp ( - ^ 

7T \ 2\k F a{p)\ 



(16) 



We probe this redefined scattering length at relative 
wavenumber p = y / 2mii/h 2 , which corresponds to the 
relative momentum of two particles on opposite sides of 
the Fermi sphere. In the regime under consideration (us- 
ing the chemical potential given in Fig. 2) the value of a(p) 
will be negative, and we find from Eq. I|I6|I basically the 
same ratio of T c /Tp as in Fig. 2. By probing the energy- 
dependent interaction at the Fermi sphere, we thus find 
in the BCS limit a result that compares excellently with 
our enhanced NSR method. 

We have studied the presently known Feshbach reso- 
nances for 6 Li and 40 K, which are the atoms currently 
used to study the BCS-BEC crossover. However, none 
of these resonances satisfy the conditions needed to ob- 
serve the results presented above. But since there are 
several bosonic systems which do satisfy the requirement 
of having a narrow resonance with large a^ gl there is no 
principal reason why a fcrmionic system could not be in 
a similar situation. Possible candidates will be fermionic 
alkali mixtures, such as 6 Li— 40 K. 

In conclusion, a positive scattering length is usually 
associated with repulsive interactions in the context of 
ultracold atomic gases. In the case of a narrow Fesh- 
bach resonance, the interactions are strongly energy de- 
pendent. We have shown that in the special case of a 
large and positive background scattering length, the in- 
teractions can be repulsive for small energies, whereas 
attractive for energies around the Fermi energy, similar 
to the effective interactions between electrons in a metal- 
lic superconductor. This results in fermionic superfluid- 
ity while the scattering length is positive, without the 
presence of bosonic molecules in the system. 

This work was supported by the Netherlands Organi- 
zation for Scientific Research (NWO). 
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